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Abstract 

Within the framework of relativistic quantum field theory, a novel method 
is established which allows to distinguish non-equilibrium states admitting 
locally a thermodynamic interpretation. The basic idea is to compare these 
states with global equilibrium states (KMS states) by means of local thermal 
observables. With the help of such observables, the states can be ordered 
into classes of increasing local thermal stability. Moreover, it is possible to 
identify states exhibiting certain specific thermal properties of interest, such 
as a definite local temperature or entropy density. The method is illustrated 
in a simple model describing the spatio-temporal evolution of a "big heat 
bang" . 

1 Introduction 



States of macroscopic systems which are in global thermal equilibrium can be dis- 
tinguished in a clearcut manner within the setting of quantum statistical mechanics 
by the Kubo-Martin-Schwinger (KMS) condition [^. This condition, which com- 
prises the pertinent features of the Gibbs ensemble in the case of finite and infinite 
quantum systems 0, has proved to be a powerful ingredient, both, in the con- 
struction of equilibrium states in concrete models and in their general structural 
analysis. 

However, there does not exist a similarly general and powerful characterization 
of states complying with the heuristic idea of being only locally in thermal equi- 
librium, e.g. having only locally a definite temperature which varies from point to 
point. The physical origin of this conceptual difficulty is obvious: if one deviates 
from a global equilibrium situation, there appear a variety of possibilities, ranging 
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from mild perturbations of equilibrium states through steady states, whose prop- 
erties are governed by external heat baths, or hydrodynamic flows up to systems 
which no longer admit a description in terms of thermodynamic notions. Thus, 
what is required are concepts allowing to distinguish between these different situ- 
ations and describing the respective local thermal properties of the corresponding 
classes of states. 

There exists an extensive literature on this subject, cf. the comprehensive treat- 
ment in P] and references quoted there. Most approaches are based on a maximal 
entropy principle, suggesting a modified Gibbs ansatz in a suitable (large time) 
limit of the theory. Other methods proceed by some procedure of "coarse grain- 
ing" directly to an effective macroscopic description. The conceptual basis of 
these approaches, however, does not always seem clear. In any case, the question 
of identifying the states of interest in the microscopic theory is left aside. 

In the present article we analyse systematically this problem in the setting 
of relativistic quantum field theory. This framework is particularly well suited 
for such an investigation since it provides an assignment of observables to space- 
time regions, yielding a description of the local properties of physical states. In 
particular, it allows a local comparison of different states, which has proved to be an 
important element of structural analysis 0. Within the present thermodynamic 
context, it was used as an ingredient in discussions of local thermal aspects of 
physical states in p, 0] and we will also make use of it here. 

Before going into technical details, let us outline the heuristic ideas underlying 
our approach. The basic ingredients are suitable "local thermal observables", 
assigned to the space-time points x, by which any given state can be compared 
with the members of the family of global equilibrium states. If a state happens 
to coincide at x with a global equihbrium state to some degree of accuracy, i.e. 
if the expectation values of a sufficiently large number of thermal observables 
at X coincide in the two states, it is meaningful to say that the given state is 
approximately in equilibrium at x. The larger the number of observables for which 
the two expectation values agree, the better the state complies with the idea of 
local equilibrium. 

The use of local observables for the description of thermal properties requires 
some justification since thermodynamic notions are usually regarded as being of 
a macroscopic nature. Let us discuss this for the case of temperature. What 
one needs to measure temperature is a suitable thermometer, i.e. a device which 
allows discrimination of equilibrium states of different temperatures. In dealing 
with very large systems in global equilibrium it does not matter where and when 
this measurement takes place. The thermometer may be coupled to any part of the 
system under consideration at any time. If one increases its sensitivity, its spatial 
extension can be taken to be arbitrarily small. Moreover, as the system is in a 
stationary state, time averages can be replaced by the mean value of sufficiently 
many measurements at arbitrarily short time intervals. Therefore, for a system in 
global equilibrium, it is possible to determine the temperature by observables in 
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arbitrarily small space-time cells, provided the corresponding measurements are 
repeated sufficiently often. The smaller the cell the larger are the fluctuations 
of the results of the individual measurements; but the collection of these data 
provides the desired information. 

This simple observation suggests the following definition of local temperature 
in a given space-time cell of a non-equilibrium system (provided a whole ensemble, 
i.e. a large number of copies of the system, can be prepared): take sufficiently 
many readings of a suitable observable in that particular cell and compare them 
with those in global equilibrium states of arbitrary temperature. If the two sets 
of data happen to coincide for some value of the temperature, it is meaningful to 
ascribe that temperature value to the chosen cell of the non-equilibrium system 
under consideration. Thus the local temperature is to be regarded as a property 
of the ensemble, not of an individual member of it. But as quantum theory is a 
statistical theory, this point of view fits naturally into its mathematical framework. 

It is conceptually convenient to proceed from the measurements in finite space- 
time cells to the idealization of temperature measurements performed at a space- 
time point. Such measurments can, however, no longer be described by observables 
(operators); one has to use quadratic forms. But since this causes no major math- 
ematical difficulties, it seems natural to perform this additional step. 

The idea of using local observables (or quadratic forms) for the description of 
the thermal properties of non-equilibrium states can be extended to other intensive 
thermal parameters of interest such as the energy density and the charge density. 
Again one has to exhibit suitable local observables which are sufficiently sensitive 
to these properties and to calibrate them by using the global equilibrium states 
as reference systems. The measured data taken in a non-equilibrium state have 
to be compared with those recorded in global equilibrium and if these data agree 
for a particular reference state one can attribute the values of the corresponding 
thermal parameters to the state describing the non-equilibrium situation. If, in 
particular, a state admits locally an interpretation in this manner in terms of a 
unique global equilibrium state, then all its local thermal parameters are fixed, 
including fundamental quantities such as the entropy density. 

Transforming this idea into action, we have to cope with a difficulty which 
we have ignored so far: the notion of equilibrium in general refers to the state 
of a system at rest with respect to the observer. Dealing with an arbitrary state 
(describing, e.g., a hydrodynamic flow) we have to take into account that these rest 
systems may vary from point to point. Consequently, the set of thermal observables 
should be sufficiently rich to determine these local rest systems and the set of 
thermal reference states should include all Lorentz-transformed equilibrium states, 
too. Moreover, in order to cover also situations, where the thermal parameters do 
not have sharp values, we have also to admit mixtures of global equilibrium states 
as thermal reference states. 

Thus, to summarize, a properly chosen set of local thermal observables together 
with the set of thermal reference states, fixed by the theory, will be our tool for 
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the analysis and interpretation of the local thermal properties of non-equilibrium 
states. 

After this exposition of the heuristic ideas underlying our approach, let us turn 
now to its mathematical formulation. In the next section we recall some basic 
elements of thermal quantum field theory and introduce the thermal reference 
states. Section 3 contains the characterization of the local thermal observables. 
In Sec. 4 we identify with their help a hierarchy of physical states which may be 
regarded as being locally close to thermal equilibrium with increasing degree of 
accuracy. An instructive model, illustrating various points of our general analysis, 
is disussed in Sec. 5 and the article concludes with remarks on some conceptual 
aspects of our approach. In the appendix we establish a theory of thermostatics 
for the set of thermal reference states used in the main text. 



2 Thermal reference states 

As we are dealing here with a conceptual problem, the present investigation is 
carried out in the general setting of thermal quantum field theory. We begin by 
recalling some relevant notions and results, primarily in order to fix our notation. 

The basic objects in our analysis are local observable fields, such as conserved 
currents, the stress energy tensor, etc., denoted generically by (with tensor 
indices omitted). Quantum fields at a point are known to be quite singular objects 
which are only defined in the sense of quadratic forms on a suitable domain, cf. the 
subsequent section. It is therefore common to proceed from the point fields 4>{x) 
to averages with suitable test functions / having compact support in Minkowski 
space M^, 

0(/) = jdxf{x)4>{x). (2.1) 

These are defined as operator-valued distributions on some common dense and 
stable domain T) in the underlying Hilbcrt space 7Y, where the theory has been 
constructed (usually the vacuum sector). The finite sums and products of the 
averaged field operators generate a *-algebra A, called the algebra of observables, 
consisting of all polynomials A of the form 

A^Y.^{fM{h)---<t^{fn). (2.2) 

We assume that also all multiples of the unit operator 1 are elements of A. As the 
fields are real (being observable), the *-operation (Hermitian conjugation) on A 
is defined by 

A*-5]0(^)...0(^)0(A), (2.3) 

where / denotes the complex conjugate of/. The Poincare transformations A G V\. 
act on A by automorphisms a\, 

ax{A) = 5^0(/i,a)0(/2,a) • • •0(/n,A), (2.4) 
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where fx is defined by f\{x) = D{\)f{\ ^x) and D is a matrix representation of 
V\_ corresponding to tlie tensor cliaracter of tlie respective field (j). If A is a pure 
translation, A = (1, a), a G M.^, we denote the corresponding automorphism by aa- 
Note that the action of these automorphisms on A is well-defined; there is no need 
to assume that they are unitarily implemented. 

As we restrict attention here to observable fields, there is one further important 
ingredient in this general setting, namely locality (Einstein causality). It is the 
requirement that fields localized at spacelike separated points commute. We will 
make use of this property at some technical points in the subsequent analysis. 

The states of the physical systems are described by positive normalized expecta- 
tion functionals on the algebra A, generically denoted by u. We recall that these 
functionals have the defining properties u^ciAi + C2A2) = ciu^Ai) + C2Uj{A2), 
ci,C2 G C (linearity), u!{A*A) > (positivity), and u^l) = 1 (normalization). 
Basic examples are states of the form 

u{A) = TipA, AeA, (2.5) 

where p is a density matrix in the underlying Hilbert space Ti and Tr denotes the 
corresponding trace. It is, however, important to note that not all states of interest 
here can be represented in this way. For, if Ti is the vacuum Hilbert space, all states 
of the form (|2.5| ) describe mixtures of local excitations of the vacuum; thus global 
thermal states are not among them. An explicit formula for the thermal states. 



analogous to ( |2.5| ), would require a change of the Hilbert space representation 
of the algebra of observables A. Since we deal here with states describing very 
different physical ensembles, it is, however, more convenient to rely on the abstract 
characterization of states given above. We mention as an aside that one can recover 
from any such state an explicit representation in terms of density matrices on some 
Hilbert space by the Gelfand-Naimark-Segal reconstruction theorem [^. But we 
make no use of this fact here. 

As has been outlined in the introduction, an important ingredient in the present 
investigation are the global equilibrium states which are characterized by the KMS 
condition [Q . Given a Lorentz frame, fixed by some positive time-like vector e G V+ 
of unit length, this condition can be stated as follows. 

KMS condition: A state satisfies the KMS condition at inverse temperature 
(3 > in the given Lorentz frame if for each pair of operators A,A'eA there is 
some function h which is analytic in the strip Sp = {z E C : < Imz < (3} and 
continuous at the boundaries such that 

h{t) = ujf3{A'ate{A)), hit + il3) = ujp{ate{A)A'), t G M. (2.6) 

In this situation ujf3 is called a KMS state. 

We will have occasion to make also use of a slightly stronger version of the KMS 
condition, proposed in . It can be established under quite natural constraints on 
the number of local degrees of freedom of the underlying theory and is a remnant 
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of the relativistic spectrum condition in the vacuum sector. It was therefore called 
relativistic KMS condition in §]. 

Relativistic KMS condition: A state U/s is said to satisfy the relativistic KMS 
condition at inverse temperature /5 > in the given Lorentz frame if for each 
pair of operators A, A' & A there is some function h which is analytic in the tube 
+ z (y+ n (/3e — V+)) such that in the sense of continuous boundary values p| 

h{a) = ujp{A'aa{A)), h{a + i(3e) = cj^(a«(A)A'), a G (2.7) 

It is a well established fact that any KMS state LVp describes an ensemble which 
is in thermal equilibrium in the distinguished Lorentz frame, describing the rest 
system of the state. As we have to keep track of both, temperatures and rest 
systems, it is convenient to combine this information into a four vector j3e G V^, 
which will again be denoted by (3. The temperature and the rest system can be 
recovered from any such temperature vector /3 G in an obvious manner. 

For given /5 G l + , the corresponding convex set C/3 of all KMS states U/3 is known 
to form a simplex under quite general conditions p|. This implies in particular 
that all states in can be distinguished with the help of "classical" (central) 
observables. Familiar examples of such observables are the chemical potential 
in theories with a conserved charge or the mean density distinguishing different 
phases. In order to simplify the subsequent discussion, we assume here that each 
family C/j is non-degenerate, i.e. for any given /3 there is a unique KMS state a;^. 
The case of more complex families C/j can be treated similarly but requires some 
heavier notation. 

Let us briefly discuss the consequences of this uniqueness assumption for the 
transformation properties of the KMS states under Poincare transformations. For 
any given KMS state cj/j and transformation A = (A, a) , the corresponding Poincare- 
transformed state is given by a;^ oa^^, where the circle denotes the composition of 
maps. It follows from the definition of KMS states that 0;^ o is again a KMS 
state belonging to the family Ca/j. In view of our uniqueness assumption we thus 
obtain the simple transformation formula 

o q;-^ = uJAfj. (2.8) 

Thus the KMS states considered here are in particular isotropic in their rest system 
and invariant under space-time translations. On the other hand, velocity trans- 
formations in general change these states, in accordance with the well known fact 
that the Lorentz group is spontaneously broken in thermal equilibrium states [^]. 
Note, however, that the temperature, fixed by the length of /3, remains unaffected 
by Poincare transformations. 

So far we have considered only KMS states, where the temperature vector (3 
is sharply defined. Yet thinking of macroscopic ensembles in thermal equilibrium 
where P is not precisely known or statistically fluctuating, it is natural to consider 
also mixtures of KMS states. We therefore define for any compact subset B dV^ 



6 



corresponding families Cb of states consisting of all mixtures of KMS states oup 
with P eB. 

In order to obtain a more explicit description of the states in Cb and to avoid 
mathematical subtleties, we assume that the KMS states cup are weakly continuous 
in i.e. all functions 



are continuous. Apart from phase transition points (which are excluded here by our 
uniqueness assumption on Cp), this property is expected to hold quite generally 
With this input, the states cub G Cb can be represented in the form 



where p is a positive normalized measure which has support in B. Note that in 
the case of degenerate families Cp the measure would also have support on the 
spectrum of all other classical observables parameterizing the KMS states. 

The elements of C = \J Cb, where the union extends over all compact sub- 
sets B C V4-, will be our reference states for the characterization of the local 
thermal observables and the interpretation of the local thermal properties of non- 
equilibrium systems. This choice turns out to be particularly convenient in the 
present investigation; yet without problems one could also admit reference states 
cub corresponding to non-compact sets B (ZV+. 

3 Local thermal observables 

Although the algebra A contains all local observables of the underlying theory, it 
is not well suited for the description of the local thermal properties of physical 
states. For, all elements of A have quite fuzzy localization properties due to 
the integration of the point fields with test functions. So they do not have an 
interpretation in terms of physical properties which can be assigned to space- 
time points. Furthermore, A contains an abundance of elements which are of no 
relevance to the thermodynamic interpretation of physical states. 

For this reason we will proceed first from the algebra A to linear spaces Qx 
of quadratic forms which can be interpreted as observables localized at the space- 
time points a; G M^. In a second step we will then select from each suitable 
subspaces T^. The elements of will be regarded as thermal observables at x 
whose thermal interpretation is fixed by the thermal reference states in C. 

3.1 Point fields as idealized local observables 

For the step from smeared field operators to point fields, let us assume that in the 
vacuum sector H there hold energy bounds of the form 



(2.9) 




(2.10) 



±0(/) < K/) (1 + i^)\ 



(3.1) 
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where / is any real test function with compact support, is some Schwartz 
norm on the space of test functions, H is the (positive) Hamiltonian, and A; > 
depends only on the field (p. This condition has a clearcut physical interpretation: 
it says that the observed values of 0(/) can become large only in states of large 
energy. Such energy bounds have been established in many models |T^. 

As has been shown in one can proceed from (|3.1| ) to an improved estimate 
(for arbitrary I > and corresponding sufficiently large m > 0) 

||(l + if)— 0(/)(l + if) — ||<c Jdx\{l-A)-^fix)l (3.2) 

where || ■ || denotes the operator norm in Ti. and A is the Laplacian on M^. It 
follows from this bound that for sufficiently large m and suitable sequences of test 
functions 5,, i G N, converging to the Dirac (5-function at the space-time point x, 
there exists the (norm) limit 

[1 + H)-"" (f){x) {1 + H)-"^ = \im{l + H)-'^(j){6i){l + H)-"'. (3.3) 

Thus the field (f){x) at x is well defined in the form sense in those states u in the 
vacuum sector, cf. (|2.5|) , which satisfy ^^((1 + HY"^) < oo. We denote the linear 
spaces of point fields (f){x), x G M^, for which ||(1 + H)~''^(f){x){l + H)^"^\\ < oo by 
Qm,x, TTL > 0. These spaces are generically finite dimensional and invariant under 
the automorphic action of the stability group of x in V]_. Evidently, Qm,x C Qm',x 
if m < m'. 

The symmetric elements of Qm,x may be regarded as (idealized) observables 
at X whose expectation values are well-defined in the states characterized above. 
Since the product of fields at a point is meaningless, their higher moments are ill- 
defined, however. We therefore propose to take as a substitute the normal products 
of these fields which can be defined with the help of the Wilson-Zimmermann 



operator-product expansion ||T2| 



A novel approach to these notions has recently been established by Bostelmann 
0. Within this general setting one finds that for each 0(x) G Qm,x and spacelike 
C G the product (f){x + ()(j){x — () is well-defined in the form sense. Moreover, for 
any given g > 0, there exist a finite number of distinguished fields (f)j{x) G Qn,x, 
n sufficiently large, and analytic functions C Cj(C), j = 1, . . . J{,(l), such that 

J{i) 

11(1 + H)--(ct>{x + 0<P{x -0-Yl ^^(0 U^)) (1 + Hr^W <c\C\' (3.4) 

i=i 

for ( tending to in any proper spacelike cone. Thus a meaningful substitute 
for the ill-defined square of 0(x) are the subspaces Af{(f)'^)q^x C Qn,x generated by 
the fields (pj{x), j = 1, . . . J{q). In a similar manner one obtains spaces M^cj)^) 



corresponding to higher powers p of 4>{x) [13|. Given r > and choosing q and n 



sufficiently large, the convergence in (|3.4| ) can be made sufficiently rapid such that 
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also 



Jii) 

1 + Hr^{d^ + c)<p{x -c)-j2d^ c,(c) M^)) (1 + I < c icr, (3.5) 



where 9^ is any monomial in the partial derivatives with respect to the components 
of (. Thus the spaces J^{(p'^)q,x contain, for large q, also fields approximating the 
"balanced" derivatives (9^ 0(x + C)(p{^ ~ C) ^oi small (. 

So far, we have restricted our attention to point fields in the vacuum sector. 
Since we are interested here primarily in thermal states, let us discuss next the 
circumstances under which these fields can be defined in such states as well. Here 
we rely on the fact that there exist also local versions of the preceding statements 
in generic cases. More precisely, in the above estimates the Hamiltonian H can 
be replaced by a suitable local Hamiltonian Hq with the help of e.g. the universal 



localizing maps considered in The operators Ho induce the time translations 
of the observables in each region C C MMn the sense that the commutators 
of Hq and of H with these observables coincide. Roughly speaking, the local 
Hamiltonians Hq may be thought of as the energy density integrated with suitable 
smoothed-out characteristic functions of compact support; so these operators can 



also be used to determine the local energy content of the states [M. Thus if a 



state u has locally finite energy in the sense that uj{{l + HoY"") < oo, we conclude 
that the expectation values of all point fields (j){x) G Qn,x, x E O, are well-defined 
in this state. 

It is evident that states allowing locally a thermal interpretation must have 
locally finite energy. So the upshot of this discussion is the insight that the states 
of interest can be extended to the spaces Qn,x, n > 0, for x varying in some region 
O. For our thermal reference states a;^ G C, this can also be seen more directly by 
noticing that, due to the invariance of these states under space-time translations, 
one has 

UBi<Pix))=MHf)) (3.6) 

for any test function / with fdxf{x) = 1. Relation (|3l6| ) may be regarded as a 
formal expression of the statement, made in the introduction, that in thermal equi- 
librium states the expectation values of intensive observables can be determined 
in arbitrarily small space-time cells. 

3.2 Selection of local thermal observables 

As has been indicated above, the algebra A, and consequently the spaces Qn,x, 
contain many elements which are of no relevance to the thermal interpretation of 
states at the space-time point x. Examples are the derivatives dx (f){x) of point fields 
which are sensitive to the spatio-temporal variations of states and therefore vanish 
in all thermal reference states. So we have to select from the spaces Qn,x suitable 
subspaces, cf. also the concluding remarks. This is accomplished as follows. 

For small n, the spaces Qn,x contain only multiples of the identity. With in- 
creasing n, there will appear, at a certain value no > in the corresponding space 
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Qno,x, also some point field (poi^) which we regard as a basic thermal observable. 
Next, as has been explained above, a meaningful substitute for the square of 4>o{x) 
are the normal product spaces Af{(f)l)q^x, q > 0. Similarly, to its higher powers 
p e N there correspond the spaces J\f{(pQ ) q,^., g > 0. Thus the elements of all these 
spaces are to be regarded as thermal observables as well. They form a proper 
subspace % of the space of all point fields, 

% = $^Ar(0o^),,., (3.7) 

p,q 

where we put J\f{(pQ) q^x = C 1 and Af{(f)Q)q^x = C0o(a;), q > 0. Assuming for sim- 
phcity that the algebra A is generated by the averages ( |2.1D of the basic observable 
00 and its normal products, we are thus led to take 7^ as the space of thermal 
observables at x. Such spaces can also be defined in more general situations in a 
systematic manner, but we refrain from giving details here. 

The above procedure introduces a natural hierarchy amongst the thermal ob- 
servables in 7^. In order to reveal this fact, let us consider any thermal ref- 
erence state ujb & C. Assuming that the KMS states satisfy the relativistic 
KMS condition, given in Sec. 2, it follows that the two-point functions ( i— » 
^^b{4>o{x =f C) 00 (a; ± C)) are the boundary values of analytic functions in the tips 
of the forward and backward tubes ± i V+, respectively. Furthermore, they 
coincide for spacelike ( because of locahty. Thus, by the Edge-of-the- Wedge The- 
orem, they are analytic for spacelike (. In view of the operator-product expansion 



( |3.4| ) and of the invariance oi lub G C under space-time translations, these two- 
point functions are therefore already determined by the expectation values of the 
observables 4>j{x) G ■Af{4>o)q,x,q > 0, in the state ub G C. In a similar manner, 
the spaces A/'(0o) q.x, g > 0, fix the p-point functions for any p E N, unless these 
functions exhibit essential singularities at coinciding points. 

As the p-point functions for large p govern the properties of those for smaller 
p (for the KMS states this is an immediate consequence of the cluster property 
which is, in turn, a consequence of their uniqueness), we conclude that the spaces 
■^{4'o)q,x, q > 0, allow, for increasing p, an increasingly finer resolution of the 
properties of the thermal states. As a matter of fact, the above argument provides 
evidence to the effect that the set Tx of all thermal observables at x in general 
separates the states in C and is sufficiently big in order to determine all properties 
of the thermal states with arbitrary precision. In other words, if one is only 
interested in the gross (macroscopic) properties of these states, it suffices to analyze 
them with the help of some subset of thermal observables taken from the spaces 
■^{4>o) q,x C Tx for small p and q. 

3.3 Macroscopic interpretation of local thermal observables 

Let us explain next how the elements of Tx provide information about the macro- 
scopic thermal properties of the states in C. In order to simplify this discussion, 
we consider only the generic case where Tx separates these states. 

As we are assuming uniqueness of the KMS states, all intensive thermal pa- 
rameters attached to these states can be represented by functions [3 t— > F{(3) of 
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the temperature vector, called thermal functions in what follows. Hence for mixed 
thermal states uob G C, the mean values of the thermal functions F are given by 



ub{F) = jdp{(3)F{(3), (3.8) 

where p is the measure appearing in the decomposition ( p.lOp oi ujb- The notation 
in ( |3.8|) indicates that the thermal functions F are to be regarded as (macroscopic) 
observables which can be evaluated in all states in C; more precisely, the states can 
uniquely be extended to these functions which appear as limits of suitable central 
sequences of local observables. 

It is crucial for our approach that the local observables in 7^ provide the same 
information about the thermal properties of the states in C as the macroscopic 
observables. Namely, one can reconstruct with their help all relevant thermal 
functions, in particular the entropy density, as is shown in the appendix. Moreover, 
since the thermal observables 7^ separate the thermal states, they can be used to 
determine the measures p in the decomposition (|2.10|) , which are needed for the 
evaluation of mean values. 

Let us exhibit this fact more explicitly. Given any (j){x) G 7^., it follows from 
the transformation formula (E.SI) that the function 



(3 ^ $(/3) = upim) (3.9) 



(which is continuous according to ( p.9|) and (|3.6| )) does not depend on x and is 



a Lorentz tensor (corresponding to the tensorial character of 0) built out of the 
temperature vector (3. It may thus be regarded as some thermal function. We use 
the fact that for each (t){x) the corresponding function $, fixing the thermal inter- 
pretation of (^(x), is known (by ( p.9|) ). As was outlined in the introduction, this 
amounts in practice to recording the mean values of the local thermal observables 
in all equilibrium states. In the following we will consistently make use of the 
above notation: lower case Greek letters, such as 0, e, denote the microscopic 
thermal observables whereas the respective upper case letters $, S, E denote the 
corresponding macroscopic thermal functions. 

Evaluating (j){x) G T^. in an arbitrary state ujb C one obtains 

uJb{<P{x)) = j dpiP) M(P{x)) = j dpiP) $(/3) = UB^, (3.10) 

which fixes some (generalized) moment of the a priori unknown measure p. By 
letting 0(x) run through all of %, the corresponding thermal functions $ run 
through a dense set of continuous functions on B C V+ (since % separates the 
states in C). Then all moments of the measure are known, hence one can recon- 
struct the measure itself from these data, cf. the subsequent section. At the same 
time this density property also implies that the space spanned by the thermal 
functions $ can be used to approximate any continuous function F on compact 
subsets B dVj^ with arbitrary precision even if F is not a member of this space. 
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This is of significance for quantities such as the entropy density, which may not 
be expected to correspond to any local thermal observable in 7^. In order to de- 
termine locally the exact mean values of such thermal functions in the states C, 
one needs information on the expectation values of an infinite number of suitable 
local thermal observables. 

On the other hand, one may expect that basic thermal functions, such as the 
thermal energy density E, can in general be determined by some element e{x) G 7^. 
Note, however, that the observable e{x) need not coincide with the full stress energy 
density 6{x) of the underlying theory even though both observables have the same 
expectation values in thermal states, 

u;p{e{x))=E{f3)=u;f,{e{x)), PeV+. (3.11) 

This peculiarity may be understood if one notices that 9{x) determines not only 
the local thermal energy of states but also other forms of energy contained, e.g., 
in the internal stress of the system. As the energy density of the states in their 
rest systems is entirely of thermal nature, the observable e{x) has to be sensitive 
only to this type of energy. This point will be exemplified in the model discussed 
in Sec. 5. 

4 Thermal properties of non-equilibrium states 

Having specified the thermal reference states and the local thermal observables 
of the theory, we have now at our disposal the necessary tools for the analysis of 
the local thermal properties of non-equilibrium states. Throughout this section we 
assume without further mentioning that the states u we are interested in can be 
extended to the local thermal observables. 

4.1 Characterization of locally thermal states 

As has been explained in the introduction, we shall determine the thermal prop- 
erties of arbitrary states by comparing them with states in C with the help of the 
local thermal observables. Given x G M'^ and any subspace <Z %;, we say that a 
state u! is iS^^-compatible with a thermal interpretation at x (jS^j-thermal, for short) 
if there exists some ujb & C such that 

u;{<P{x)) = cub(0(x)), 0(x) G S,. (4.1) 

Thus u looks like a thermal state with respect to all observables in Sx. Under 
these circumstances one can consistently define mean values of the corresponding 
thermal functions in the state u at the space-time point x, setting 

o;(<l>)(x) =^(0(a;)), 0(a;) G 5,, (4.2) 

where $ corresponds to 0(x), cf. (|3l9|). For if two observables (f)i{x) , (f)2{x) G % 
give rise to the same $, it follows that ujB{(t)i{x)) = ujb{^) = ^^B{(t>2{x)) for all 
ujb G C. So, for the given x, the basic relation ( [4.2| ) provides a lift of uj to the 
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subspace of thermal functions $ fixed by S^- It thereby leads to a local thermal 
interpretation of this state. 

The preceding characterization of states admitting locally a thermal interpreta- 
tion is physically meaningful but difficult to use in practice because of the apparent 
need to compare these states with all members of the family C. It is therefore grat- 
ifying that there exists an equivalent characterization which relies entirely on the 
space of thermal observables To see this we introduce on this space a family 
of seminorms tb for all compact subsets B C V+, setting 

rB(0(x)) =sup mf3)\. (4.3) 

Bearing in mind that the functions $ are continuous, it is clear that the supremum 
exists. The following criterion makes use of this notion. 

Criterion: Let Sx be any subspace of 7^ containing the identity and let u be any 
state on A. The state is iS^^-thermal if and only if there is some compact subset 
B C V+ such that 

<rs(0(a;)), <j){x) e S,. (4.4) 



The latter condition has a simple physical interpretation: the mean values of 
the local thermal observables should not exceed their maximal possible values in 
the thermal states. This constraint can be checked more easily in applications. 

The mathematical justification of the criterion relies on standard measure theo- 
retic arguments: let u be iSa.-thermal. Then there exists some positive, normalized 
measure p with support in B, cf. ( |2.1CI| ), such that 



u;icf>{x)) = coB{(pix)) = / dpiP) co^i^ix)) = / dp{P) $(/?), 0(x) G S,. (4.5) 



Relation ( f4.4| ) then follows by a straightforward estimate. Conversely, if relation 
( [4.4|) holds for u and some B, the linear functional uj{^){x) = uj{(f){x)) on the 
space of thermal functions $ corresponding to 0(a;) G Sx is well-defined. For $ = 
implies TB{(t){x)) = sup^g^ l'^*(/5)l = and consequently a;($)(x) = 0. Moreover, 
as |co'(<I>)(a;)| < tb{(P{x)) = sup^g^ l^(/5)|! this functional can be extended in $ by 
the Hahn-Banach Theorem to the space of all continuous functions F on B such 
that 

\uj{F){x)\<snp \F{(3)\ = \\F\\b. (4.6) 

l3eB 

As the functions F with norm 1 1 ■ 1 1^ form a C*-algebra with unit and as uj{l){x) = 
1, it follows that uj{ ■ ){x) defines a positive, normalized functional on this algebra 



[0. But any such functional can be represented in the form |]T6 



coiF)ix)= dpiP)FiP), (4.7) 
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where p is some positive, normalized measure with support in B. Setting 



ujb{A) = Jdp{(3)ujp{A), AeA, {U 



we have thus found a state a;^ G C which coincides with u in the sense of condition 
Llf). Hence u is iSr-thermal. 



The thermal states ub describing the local thermal properties of u are not fixed 
by condition (|4.1j ) if the set of thermal observables Sx is too small. If u admits, 
however, a thermal interpretation in the sense of ([4.1| ) with respect to all thermal 
observables %, then the corresponding state ub & C (and hence the corresponding 
measure p) is in general unique by the arguments given in the preceding section. 

4.2 Existence of locally thermal non-equilibrium states 

For any given finite dimensional subspace Sx of local thermal observables and 
any compact subset B <Z V+, there exists an abundance of states u on A which 
coincide on with some thermal state ub G Cb but are not 7^-compatible with 
a thermal interpretation. So they correspond to a non-equilibrium situation at 
X which admits, however, an interpretation in terms of the subset of thermal 
functions corresponding to Sx- 

To establish this fact, let us first assume that tb defines a norm on Sx, i.e. 
TB{(f>{x)) = for (f){x) e Sx implies 0(x) = 0. Picking any state uq on A and 
taking into account that on a finite dimensional space all linear functionals are 
continuous and all norms are equivalent, we get 

M(^{x))\ < ctb{(P{x)), 0(x) e Sx (4.9) 

for some suitable constant. Thus we can lift ujq to the space of thermal functions, 
setting a;o($)(x) = Uo{(f){x)), (f){x) G Sx. Moreover, by the Hahn-Banach Theorem, 
the latter functional can be extended to the space of all continuous functions F 
on B such that 

\uJo{F){x)\<c sup \F{f3)\. (4.10) 

We may assume that this extension is hermitian, replacing ujq{F){x) by the ex- 
pression 2~^(ti;o(-F)(x) + u!o{F*){x)), if necessary. As any linear, hermitian and 
continuous functional on the C*-algebra of continuous functions on B can be rep- 
resented by a signed measure a with support in B [|16| , we conclude that 



uJo{F){x) = jda{(3)F{(3). (4.11) 
Decomposing a into its positive and negative parts a-t, o" = (T+ — (t_, and setting 
u{A)={l + a_{B))-\uj^{A) + Jda4P)ujf,iA)), A E A, (4.12) 
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we have thus exhibited a state uj on A such that 

cu(0(x)) =cub(0(x)) = a+(fi)-i [da+iP)cu^i(Pix)), 0(x) G 5,, (4.13) 



i.e. UJ is iSa;-thermaL But in general it is not 7^-thermal since c^o was completely 
arbitrary. As is seen from relation ( [4.12|) , the state u so constructed may be 



interpreted as a perturbation of a thermal background state in Cb- 

Let us turn now to the case where there is some non-trivial subspace Sx^o C Sx 
which is annihilated by tb- If, as expected in the absence of phase transitions, 
the thermal functions $ are analytic on V+, the relation TB{(f){x)) = implies 
$ = if -B has an open interior, i.e. there exist non-trivial relations between 
the thermal functions $, (f){x) G Sx- The existence of such relations has a simple 
physical interpretation: it amounts to the existence of equations of state (such as 
the relation between energy density and pressure, for example). 

For the construction of non-equilibrium states which are jS^^-thermal, we pick 
now any state ujq on A which complies with these equations, i.e. which annihilates 
the subspace Sx^q. Such states which are not 7^-thermal also ought to exist in 
abundance; for the local validity of some equation of state does in general not 
imply that the system is locally in perfect thermal equilibrium in the sense that 
all correlations described by the elements of 7^ are of a thermal nature. We will 
exemplify this fact in the model discussed in the subsequent section. 

In view of its proper kernel, the functional Uo can be projected to the (finite 
dimensional) quotient space Sx/Sx,o on which tb induces a norm. Hence, as be- 
fore, Uq satisfies the bound ( |4.9| ) and the subsequent construction can be carried 
out. So we conclude that there exist non-equilibrium states admitting a thermal 
interpretation for any given finite dimensional subspace of thermal observables. 

4.3 The degree of thermal stability 

Let us discuss next how the states admitting locally a thermal interpretation can 
be ordered into classes of increasing thermal stability (thermalization). Here enters 
the hierarchical structure of the thermal observables 7^, exhibited in the preceding 
section. Starting with the trivial subspace Sx = CI, all states are Sx-thermal since 
only their normalization is tested. By adding to Sx the basic thermal observable 
4>o{x) one imposes already some non-trivial constraint. For (j)o{x) is an unbounded 
quadratic form, hence arbitrary states generically violate condition ( |4.4| ) for given 
B. Next, one adds to Sx the elements of the normal product spaces J\f{4>l) q,x for in- 
creasing q which provide information about the correlations of the basic observable 
4>o{x) at neighboring points. Using these enlarged spaces Sx in the compatibility 
condition (^4.1|) one selects subsets of states where these correlations are of a ther- 
mal nature. Note that these constraints resemble the conditions imposed on states 
in the derivation of transport equations. 

Proceeding in this manner, one arrives at spaces Sx containing also elements 
of J^{(f)o) q^x for higher powers p. As has been explained in the preceding section, 
the resulting compatibility conditions distinguish states which exhibit increasingly 
more refined features of thermal equilibrium states and, in this sense, come closer 
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to the idea of a genuine equihbrium situation. Thus, for states which are S^- 
thermal, the size of may be taken as a measure for their degree of thermal 
stability. 

4.4 Determination of specific thermal properties 

To judge whether in a given state uj some thermal function $ has locally a def- 
inite value, one has to make sure that the state is compatible with a thermal 
interpretation on some sufficiently large space S^- For one has to determine not 
only the mean value of $ but also its statistical fluctuations. If, for example, 
contains observables 0i(x) and 02 (a^) corresponding to $ and respectively, then 
the observable 

5(f)^{x) = (f)2{x) - 2k (j)i{x) + 1, K G M, (4.14) 

corresponds to the thermal function ($ — /t 1)^. It is therefore non-negative in all 
thermal reference states C and vanishes only in those states in which $ has the 
sharp value k. Thus if u is iSa;-thermal and uj{6(j)K{x)) = for some k, we may 
conclude that $ has the sharp value At at a; in this state. It is apparent that this 
conclusion can also be drawn under more general conditions. 

For suitable spaces Sx one can also distinguish in this manner states which 
locally have a sharp temperature vector j3, i.e. which coincide on Sx with the KMS 
state Up. The minimal spaces admitting such an analysis are finite dimensional in 
generic cases, cf. the discussion in the subsequent section. 

If a state u coincides on Sx with some KMS state ujp, all thermal functions $ 
corresponding to 0(x) G Sx have locally definite values in this state. However, by 
enlarging Sx, w may cease to have a thermal interpretation. Phrased differently, 
UJ may share only certain gross thermal properties with the KMS state and a 
more refined analysis would reveal its non-equilibrium nature. Thus the method 
of analyzing the local thermal properties of states with the help of suitable subsets 
of thermal observables amounts to some procedure of coarse graining. It leads to 
the identification of an abundance of states having locally certain definite thermal 
properties. Such an approach seems natural if one deals with non-equilibrium 
systems. 

4.5 Space-time evolution of thermal properties 

The formalism established so far can easily be extended to states admitting a 
thermal interpretation in subregions (9 C M^. Yet in order to obtain a sufficiently 
simple description, one has to keep the thermal functions fixed in the respective 
regions. This is accomplished by identifying the spaces Sx, x & O, with the help 
of the automorphic action of the translations, setting 

Sx = axiSo), xeO. (4.15) 

With this convention understood, we say a state u is iSo-compatible with a thermal 
interpretation in O if for each x G O there is some tu^ G C (in general depending 
on x) such that relation ( [4.1| ) holds for the respective Sx. The resulting functions 

x^u{<t>){x) =uj{(f){x)), xeO, (4.16) 



16 



then describe the space-time behavior of the mean values of the thermal functions 
$. Hence they provide the link between the microscopic dynamics and the evolu- 
tion of the macroscopic thermal properties, i. e. the thermodynamics of the states. 
It is an intriguing question under which circumstances this evolution can also be 
described in terms of transport equations. 

We have thus solved the conceptual problem of identifying non-equilibrium 
states admitting locally a thermal interpretation, and of describing their specific 
thermodynamic properties. This formalism can now be applied to the analysis 
of non-equilibrium states in concrete models, as will be exemplified in Sec. 5. 
Furthermore, it is a suitable starting point for a general structural analysis of 
these states. 

4.6 Field equations and quasiparticles 

As a first illustration of the latter point, let us show how the microscopic dynamics 
leads to linear differential equations for the evolution of mean values of thermal 
functions in states which are sufficiently close to thermal equilibrium. We start 
from the assumption that our basic observable (j)o{x) satisfies a field equation of 
the form \I\x(f)o{x) = C,o{x), where C,o{x) is some linear combination of elements in 
the normal product spaces J\f{4>Q°) qo,x for certain pq and go (hence it is a thermal 
observable). As uJi3{^o{x)) = upiDxCpo^x)) = DxUJi3{(I)q{x)) = because of the 
invariance of the KMS-states under space-time translations, we conclude that the 
thermal function Hq corresponding to ^oi^) vanishes. 

Now let Sx be any space of thermal observables containing the above normal 
product spaces J^{4>o'^) go,x, hence in particular 0o(x) and ^o{x), and let cu be any 
state which is iSo-compatible with a thermal interpretation in O. Then the local 
mean values of the thermal function $o in this state necessarily satisfy the wave 
equation in O, as is seen from the chain of equalities 

nxuj{%){x) = DxUjiMx)) = uj{nxMx)) = ^(eo(a;)) = uj{Eo){x) = 0. (4.17) 



There exist many other thermal observables for which such behavior of the 
mean values in suitable states can be established. To reveal the underlying simple 
mechanism, we proceed from the elementary relation 

= -De Mx + C)0o(x - C) + 2eo(x + C)0o(x - C) + 20o(x + CM^ - C), 

which holds in the form sense for spacelike (. Taking into account the remark after 
relation ( |3.5|) about balanced derivatives, it is clear that the expressions on the 
right hand side of this equality can be approximated in the limit of small spacelike ( 
by elements of the normal product spaces M^cj)^ q^x for suitable p and q. Similarly, 
the expression (pQ^x + Q't'oi.^ ~ C) can be approximated by elements of Af{(j)l)q^x- 
Comparing both sides of relation ( ^.18| ), one thus finds thermal observables (/'(x) 



in J^{(pQ)q^x for which \Dx(j){x) = ^(x) is also a thermal observable. As before one 
can then show that the local mean values of the thermal function $ corresponding 
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to 0(a;) are solutions of the wave equation in states which are compatible with a 
thermal interpretation on sufficiently large subspaces of thermal observables. 

So the space-time evolution of these thermal functions exhibits patterns of a 
massless particle propagating through the state a;, provided oj is locally sufficiently 
close to thermal equilibrium. Note that this thermal compatibility condition im- 
poses quite stringent constraints on the state in the case of interaction, since the 
spaces Sx then have to contain thermal observables in M'{4>o)q,x for p > 2. For 
states shghtly violating this condition in the sense that u;{^{x)) is different from 
but small, the above equations are, however, still valid in an approximate sense. 
This result is in accordance with the familiar quasi-particle interpretation of per- 
turbations of equilibrium states. It emerges here as a by-product of our approach 
to the characterization of the local thermal properties of non-equilibrium states. 

5 An instructive example 

In this section we illustrate the preceding abstract notions and results in the theory 
of a free massless scalar field. We have chosen this particularly simple example 
since it allows the elementary computation of many quantities of interest by scaling 
arguments. After a brief outline of the model we will determine the structure and 
physical significance of its local thermal observables. Wc then exhibit interesting 
examples of non-equilibrium states, describing a "big heat bang", for which a 
definite temperature, thermal energy and entropy density can be defined at every 
space-time point in the future cone of some given initial point. 

5.1 The model 

The free massless scalar field (f)o{x) on R^, playing the role of a basic thermal 
observable in the present model, is characterized by the field equation and com- 
mutation relation 



n,Mx) = 0, [Mxi):Mx2)] = (27r)-3 / rfpe-^(-i--^)f£(po)5(/) ■ 1- (5.1) 



It generates a polynomial *-algebra A, describing the local observables of the 
theory. This algebra is stable under the actions of the Poincare group V]., given 
by a\^a{4>o{x)) = (poi-^x + ^t)) the dilations M+, given by Ss{(f)o{x)) — s(f)o{sx), and 
the gauge group Z^, given by j{(f)Q{x)) = —(()q{x). 

We restrict attention here to states uj on A which are gauge invariant, i.e. 
(^07 = cu, so the respective n-point functions of (po vanish if n is odd. The simplest 
examples of this type are quasifree states. They are completely determined by their 
two-point functions through the formula 



LLl{4>o{Xi)(j)Q{x2) ■ ■ ■ (po{Xn)) 

^ /EpaiH„g.'^('/'o(a;n)0o(5 





)Mxin)) ^ even. 



(5.2) 



n odd. 
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Whenever a (generalized) function u;(0o(a;i) 00(2^2)) is consistent with the con- 
straints imposed by (|5.1| ) and the positivity condition i-t;(0o(/)0o(/)) > for all 
test functions /, the functional uo obtained by linear extension from ( p. 21 ) defines 
a state on A. Thus there is a large supply of simple states from which we will also 
draw the non-equilibrium states considered below. 

5.2 Thermal reference states 

It is a well known fact that the algebra A has a unique gauge invariant KMS state 
Up for each temperature vector [3 eV^. This state is quasifree, so it is determined 
by its two-point function given by 

u;p{U^,)U^,)) = (27r)-3y"rfpe-'(^-^^)%(po)5(p')Y3^. (5.3) 

We mention as an aside that these states satisfy the relativistic KMS condition 
and comply with our continuity assumption ( |2.9| ). In fact, the respective functions 
are analytic on V^, as is easily verified. 

As outhned in Sec. 2, the KMS states up, /3 G V+, fix the convex set C of 
thermal reference states which enters into our analysis of the thermal properties 
of (gauge invariant) non-equilibrium states. 

5.3 Local thermal observables 

Let us turn now to the analysis of the local thermal observables in this model. 
We will consider primarily the spaces A/'(0o) g,x, Q > 0, which are generated by the 
Wick square of (f)Q, its balanced derivatives and the unit operator 1. Introducing 
the multi-index notation = (/ii, fi2, ■ ■ ■ /^m) and setting = d^^^ ' ' ' '^(^ ' these 
balanced derivatives are defined by 

: 02 : (a;) = lim {M^ + QM^ - C) " ^00 {M^ + OM^ - 0) l) , (5-4) 

where cUqo denotes the vacuum state (which can be recovered from (|5.3|) in the 
limit of large time-like jS). Note that for odd m the balanced derivatives vanish, 
since (po^x + C)4>o{^ ~ C) is even in ^ as a consequence of locality. 

For the determination of the thermal functions corresponding to these observ- 
ables one proceeds from the relation 

^/siM^ + OM^ - 0) - t^oo(0o(a; + 000(2; - c)) 

= (2vr)-^|^ cos(2C|>) ie^^^-ir' = (2^)"'E/ p cos(2Cp) e^^E , ^^'^^ 

where p = {\p\,p)- The even derivatives with respect to (, when applied to the 
latter integral, can be replaced by d^, multiplied with appropriate constants. We 
omit the remaining simple computations and only state the final result: 
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where = if m is odd and = (— l)™'/^(47r)™(m + 2)\^^Bm+2 if is even, 
Bn being the Bernoulh numbers. It is instructive to have a closer look at these 
functions for small m. 

(a) The thermal function attached to the density :(^q: {x) is jS ^ (12/5^)"^, i.e. the 
square of the temperature (apart from a constant). 

(b) The balanced derivatives e'^'^(x) = —(1/4) Q^^" -.(pQ-. (x) give rise to the thermal 
functions 

P ^ E'^'iP) = (ttVoO) {AP^'P" - /? V") (/?')"^- (5.7) 

So they coincide with the expectation values of the components of the standard 
(symmetric and traceless) stress energy tensor 9^^{x) in the states a;^. This fact 
can be understood if one notices that this tensor can be represented in the form 

e>^\x) = e^^(x) + (1/12) {d^,d: -g^^n,) :0^(x), (5.8) 

where the second tensor on the right hand side vanishes in all equilibrium states 
because of the derivatives with respect to x. Since this tensor is the coboundary of 
the vector-valued two-form (1/12) (9^ g'"-' — g^^) : 0q : (x), it follows from Gauss' 
law that it also does not contribute to the total energy of states which deviate 
only locally from an equilibrium situation. On the other hand, it is of relevance 
in non-equilibrium states, where it describes the transport of energy driven by 
sources which are localized at the boundary of the system. 

We call e^^{x) the thermal energy tensor. It is conserved and symmetric, but 
its trace e^(x) = (1/4) Da. : 0q : (x) does not vanish identically. According to the 
discussion at the end of the preceding section, it is, however, zero in all states 
which are sufficiently close to equilibrium. 

The thermal energy tensor and the Wick square can be used to distinguish KMS 
states in C corresponding to a given temperature vector. To establish this fact we 
consider, for fixed x G V"+, the thermal observable 

5^e{x) = (30/7r2)x^x,e'^"(a;) - 24 :0^(a;) + (x^)-! 1 (5.9) 

with corresponding thermal function 

P ^ (l/3)(4(/3x)2 - p^x'){P''Y^ - 2{P^)-^ + (x2)~i. (5.10) 

As {P>i)'^ > where equality holds if and only if /5, x G V+ are parallel, this 

function is evidently positive on V+ apart from the point /3 = x, where it vanishes. 
Hence if UB{S^e{x)) = for some state cus G C one can conclude that it must be 
the KMS state cu^^ corresponding to the temperature vector x. Thus in this model 
a finite number of thermal observables are sufficient to decide whether a state has 
locally a definite temperature and a well-defined rest system. 

(c) Turning to the higher balanced derivatives, it follows from (|5.6| ) that all corre- 
sponding thermal functions are solutions of the wave equation since {P"^)^^ = 
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on V+. So the subspace generated by these functions does not separate the states 
in C. Yet it is sufficiently large in order to approximate the entropy current S^^, 

f] ^ S^'{(3) = (27rV45) (5.11) 

on any compact subset B G V^. (The determination of the entropy current from 
the microscopic data is outlined in the appendix.) 

For the proof of this assertion we pick any lightlike vector / and contract the 
balanced derivative : 0o ■ (^) with the tensor l^. The corresponding thermal 
function is, cf. relation (|5.6| ), 

/? ^ c„ {W^r (/?')"' = Cm m\ 2"^ {lf3r{(3'r"'-' (5.12) 

if m is even; for odd m it is identically zero. As {Ij3){j3'^)~^ is positive and bounded 
on any compact subset B dV^ and the (positive) square root is holomorphic on 

the right complex half plane, it follows that /3 ^ (//3)(/3^)"^ = {{WfiP^y^y^ 
can be represented on i? as a convergent power series involving only even powers 
(Z/3)™'(/5^)~'", m G 2No. This proves that (3 y-* IS{(3) can be approximated on B by 
linear combinations of the thermal functions (3 i— > (//5)'"(/3^)~''"~^, m G 2No. Since 
/ was arbitrary and the lightlike vectors generate a basis of R^, we conclude that 
also the entropy current can be approximated on any B with arbitrary precision by 
the thermal functions corresponding to the balanced derivatives. Thus any state 
which admits a thermal interpretation on the space generated by these derivatives 
has a well-defined mean entropy density at the space-time point x. 

5.4 Examples of non-equilibrium states 

We have already seen in the abstract analysis that there exist non-equilibrium 
states admitting a thermal interpretation on finite dimensional subspaces of ther- 
mal observables. Here we consider non-equilibrium states which are substantially 
closer to a thermal equilibrium situation, namely, these states admit a thermal 
interpretation on the infinite dimensional spaces Sx of thermal observables gener- 
ated by A/'(0o) q^x, p = 0, . . . 3, g > 0, for all x in some future hght cone. Moreover, 
they have a definite temperature, thermal energy and entropy density at all of 
these space-time points. 

We are taking here the simplest states of this type in order to illustrate our 
general method. Yet these examples are also of some physical interest. They 
describe the spatio-temporal evolution of systems which have infinite temperature 
at some space-time point (corresponding to a "big heat bang"). Although we 
are dealing with a massless free field theory, our results provide some idea of the 
dynamical effects of such singularities in more realistic theories. For the masses of 
particles and their interaction should play a secondary role in the neighborhood 
of such singular points. 

The non-equilibrium states considered here are quasi-free with two-point func- 
tions of the form, 7 > being fixed and Xi, X2 G l + , 

ujUMxi)Mx2)) = (2vr)-^y"rfpe-(--^-^)%(po)5(p')Y3^z4;7T^;)^- (5-13) 
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These functions are consistent with the field equation and commutation relation 
in They are also invariant under Lorentz transformations and dilations but 

not under translations. Moreover, 

t^bhb(0o(a;i)0o(a;2)) 

(5.14) 



[ dp 



hence, decomposing the cosine according to Euler's formula, it follows at once 
that each functional cUbhb satisfies the positivity condition ci;bhb(0o(/)0o(/)) > if 
supp/ C V+. By standard arguments it can therefore be extended to all of ^ to a 
(singular) state. Since we are not interested here in this extension we do not dwell 
upon this point any further. 

The computation of the expectation values of the balanced derivatives yields 
for even m and x G VO. 



a;bbb(S'':0^(x)) 

= (27r)-3f; /■^(_i)W22™pMe-"2^^^ = ^,(.) (5/^ : 0^ (x)) , 

n=l 



(5.15) 



where oj^^x) is the KMS state corresponding to the temperature vector fi{x) = 2'jx. 
Thus ciJbhb is -thermal, Sx being the subspaces of thermal observables generated 
by A/'(0Q)g,a;, p = 0, ...3, g > 0. Note that tUbhb is compatible with a thermal 
interpretation on the spaces A/'((/>o) q^x for odd p and g > since it is, hke the KMS 
states, gauge invariant. It has, however, no longer a thermal interpretation on 
the spaces Af{(pQ)q^x, g > 0, as one verifies by direct computation. So the higher 
correlations of the field 0o in this state are of a non-thermal nature. But the state 
approaches equilibrium for large time-like translations a in the sense that for all 
local observables A ^ A 

lim tUbhb o aa{A) = Uoo{A), (5.16) 

a 

i.e. ciJbhb looks asymptotically like the vacuum cUoo- 

Let us turn now to the local thermal interpretation of this state. As 5j^e(2:) e S^, 
cf. relation (|5.9|), and uj^i,^[5^e{x)) = for x = P{x), we conclude that cUbhb 
coincides on Sx with the KMS state x E V+ (which is also apparent from 

relation ( |5.15| )). Its local temperature T{x) = (/3(x)^)~^''^ = (47^x^)~^/^ decreases 
monotonically in time-like directions and all Lorentz observers moving along some 
world line M+ e register the same temperature at a given time after the big heat 
bang at 0. 

The thermal energy density in the state a;bhb at the space-time points x E V-^ 
can be read off from the tensor 

^b..b(E'^'^)(x) = uj,,,{e''''ix)) = iny9Q){AP''{xW{x) - g^''Pixf)il3ix)Y' 

= (7rVl4407^)(4a;^x'^ -^'^V)(x2)-^ (5.17) 
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describing a flow of massless particles in which is isotropic for the above Lorentz 
observers. They also find that the relation between the thermal energy density 
and temperature is in accordance with Stefan-Boltzmann's law for scalar massless 
particles. 

Since e^(x) = (1/4) : 0o • (^) because a;bhb(E^)(a^) = 0, it follows that 
ciJbhb(: 00 • (^)) — 0' so the density : (j)"^ : (x) propagates through the state like 
a massless particle, in accordance with the general results in Sec. 4. As a matter 
of fact, relation ( |5.15| ) implies that D^,. tUbhb (S'^ '■ 4>o '■ {^)) = for all balanced 
derivatives. 

It is an interesting fact that the full energy density of the state cUb^b is larger 
than its thermal energy density. Making use of relation ( |5.8| ) one obtains 

uJ^U^^'^ix)) = ((7rVl4407^) + (I/28872)) (4x^x^ - g^^'x^) {x^^. (5.18) 

Thus the coboundary term in relation ( ^.8|) leads to an additional contribution 
due to the transport of energy from the hot boundary of the light cone into its 
interior. Note that for 7 = 1 this term is of the same order of magnitude as the 
thermal energy, but it is not visible in the local thermal properties of the state. 
In other words, if one would try to determine the energy density of the state from 
its temperature by using Stefan-Boltzmann's law, one would underestimate its 
value. This feature of an apparently "missing energy" may be of significance in 
cosmological models. 

As '-^bhb is compatible with a thermal interpretation on all spaces Af{(f)Q) q,xi 
q > 0, one can consistently attribute to it the mean entropy current 

uJ^iUS^'){x) = (27rV45) (3^'{x){(3{xf)-^ = (iryiSO-f^) x^'ix^)-^ (5.19) 

for X E cf. the discussion in the preceding subsection. Since the model does 
not describe dissipative effects due to the lack of interaction, it does not come 
unexpectedly that this current is conserved, 

d,uj,,,{Sn{x) = 0, (5.20) 

so there is no entropy production. The entropy current decreases monotonically 
in time-like directions. But the total entropy within spacelike sections of of 
the form \x\ < v\xo\ for fixed f < 1 stays constant for xq > 0. Thus, from a 
macroscopic point of view, the state cUbhb describes an equilibrium situation at all 
points in V+ in spite of the fact that it is microscopically out of equilibrium. 

We conclude our discussion of the state cjbhb with the remark that one can 
generate from it other states which are still 5y+-thermal. Namely, for any given 
translation a G V+ the state cUbhb o «a has this property, too, since aa{Sv+) C Sy^ 
and the thermal reference states in C are invariant under translations. Moreover, 
since C is stable under convex combinations, the states / du^a) cjbhb ° «a, where 
1/ is any positive normalized measure with compact support in V+, are also Sy^- 
thermal. If u differs from the Dirac measure, these states have, however, locally 
no longer a definite temperature vector and the corresponding expectation values 
of the thermal functions exhibit a more complex space-time behavior. 
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6 Concluding remarks 



We have established in the present article a novel method for the characterization 
and analysis of non-equilibrium states admitting locally a thermal interpretation. 
The basic idea in our approach is to compare these states with global equilibrium 
states by means of local observables. The inevitable step of "coarse graining" 
in the passage from a microscopic theory to a macroscopic description has been 
accomplished by restricting attention to distinguished subspaces of thermal ob- 
servables, in accordance with the basic ideas of statistical mechanics having their 
origin with Boltzmann. As the selection of these observables is of vital importance 
in our approach, it seems appropriate to comment on this point in somewhat more 
detail. 

At first sight it might seem natural to implement the idea of comparing states 
by defining a suitable distance between them. This can be accomplished by intro- 
ducing a norm || ■ || on the observables A {e.g. by fixing suitable Schwartz norms 
on the underlying spaces of test functions) and making use of the inherent net 
structure of the observables, providing to each space-time region (9 C some 
subalgebra A{0) C A generated by observables A whose underlying test functions 
have support in O, cf. relation ( ^.21 ). A physically meaningful local norm distance 
between any two states u!i,u!2 can then be defined, setting 

||cji -CJ2II0 = snp {\iUi{A) - uj2iA)\ : A e AiO), \\A\\ < 1}, (6.1) 

which is just the maximal possible difference of the expectation values of normal- 
ized observables in A{0) in the two states. 

Making use of this notion, one could try to distinguish non-equilibrium states 
uj admitting a thermal interpretation at space-time point x by the condition that 
there is some state ujb G C and some g > such that for each sequence of regions 
Og 3 X with diameter e 

\\uJ -ujbWo, < ce'^ for e ^ 0. (6.2) 

For states u coming close to the idea of local equilibrium, the exponent q should 
be large. This condition turns out to be too stringent, however. Namely if uj has 
locally finite energy, it follows from ( |6.2| ) that u and ub coincide on some space of 
quadratic forms Qn(q),x, where n{q) is large for large q. But for sufficiently large n 
the spaces Qn,x contain, together with some observable 4>{x), also its derivatives 
dx4>{x) |jl3[- Hence dxUj{(j){x)) = uj{dx(t>{x)) = u!B{dx(p{x)) = if u is, in the 
above sense, sufficiently close to local thermal equilibrium. So, for states satisfying 
the above condition in some space-time region, the expectation values of these 
observables necessarily exhibit a trivial spatio-temporal behavior in that region. 
In other words, non-equilibrium states which are only locally close to equilibrium 
in general do not comply with this strong constraint. 

It is noteworthy that one arrives at a similar conclusion in the C*-algebraic 
setting of local quantum physics, where the forms Qn,x appear as dual spaces of 



the germs of states of (locally) finite energy |T3|, |T 
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We have solved this conceptual problem by distinguishing certain subspaces 7^ 
of thermal observables, relying on the notion of normal product. Intuitively, these 
observables describe density-hke quantities which can be attributed to space-time 
points, in contrast to their derivatives, which are sensitive to the spatio-temporal 
changes of states. The spaces 7^ have a natural hierarchical structure, ranging 
from observables which are sensitive to gross thermal properties up to quantities 
which allow to determine the subtle features of higher correlations. 

Non-equilibrium states admitting a thermal interpretation can be distinguished 
by the condition that they coincide with thermal reference states on subspaces of 
the thermal observables and their degree of thermal stability can be read off 
from the size of the respective spaces. This characterization is the appropriate 
substitute for the apparently too restrictive condition ( |6.2|) . 

As has been explained in Sec. 4, the spaces still contain certain derivatives 
of observables, but their appearance is of dynamical origin (it is a consequence of 
the field equations). Their existence leads, on the one hand, to relations between 
the thermal functions, corresponding to equations of state. On the other hand, 
they explain the quasi-particle-like propagation of densities in states which are 
sufficiently close to thermal equilibrium. 

We have restricted our attention to the simplest class of theories, both, with 
regard to the structure of the observables and the family of thermal reference states. 
But our arguments can be extended to more complex situations with only little 
more effort. We therefore believe that our approach provides a natural setting for 
the analysis of the thermodynamic properties of non-equilibrium states in concrete 
models as well as in the general framework of quantum field theory. It should help 
to shed new light on the complex features of non-equilibrium systems. 

A particularly interesting issue is the phenomenon of thermalization (approach 
to equilibrium) which is related to the problem of the arrow of time. As any local 
quantum field theory is PCT-symmetric and the spaces of thermal observables 
transform covariantly under this symmetry, there is no such arrow encoded in the 
algebraic structures used in the present approach to select the states admitting 
a thermal interpretation. Yet the states can break the PCT symmetry in the 
sense that they are compatible with a thermal interpretation either in some future 
or in some past light cone, but not in all of Minkowski space (unless they are 
in global equilibrium). The state cUbhb considered in Sec. 5 nicely illustrates this 
point: it does not have a thermal interpretation extending beyond V+, and similarly 
the state uj^^^ o t?, obtained from tUbhb by the action of the anti-automorphism 
implementing the PCT symmetry, is thermal only in the past light cone — V^+. A 
proof establishing such a one-sidedness of thermalization in generic cases would 
be an important step towards the understanding of the microscopic origin of the 
arrow of time. We hope to return to this problem elsewhere. 
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A Thermostatics in C 



In this appendix we determine the thermostatic properties of the thermal reference 
states. As in the main text, we assume that for each temperature vector j3 G 
V+ there is a unique KMS state ujp. By apphcation of the fundamental laws of 
thermodynamics, the thermal functions can be determined if the corresponding 
expectation values of the stress energy tensor 9^^{x) or of the thermal energy 
tensor e^'^(x) G 7^, respectively, are known [|T^. We will discuss in which sense 
these functions retain their interpretation in thermal reference states which are no 
proper equilibrium states. 



In view of the tensor character of 6^'^{x) and relation (2^), the corresponding 
thermal function E'^'^ has the form 

(3 ^ = uj^ie^^ix)) = Q{f3^) e^e'' - P{(3^) g^"", (A.l) 

where e = specifies the rest system of the state and Q, P depend on the 

underlying theory. As we are dealing with an equilibrium situation, the function 
P can be interpreted as pressure and Q{P^) = — 2/5^P'(/3^) [jT8|, Ch. 4]. Because of 



relativistic covariance, the entropy current is of the form 

f3^ S^'ip) = S{(3^)ef', (A.2) 



where S is fixed by the Gibbs relation in the rest system ||T8|, Ch. 4], 

SW') = WY'QiP')- (A.3) 

The free energy density F'^'^ is given by the tensor 

P ^ F^'-'iP) = E^'ip) - {l/2){PY'^\e''S''iP) + e'S^iP)). (A.4) 

Inserting the preceding expressions for the entropy and energy, one finds that it 
coincides in all Lorentz frames with the negative pressure, 

F'"'{(3) = -P{p^) g^'r (A.5) 

Turning to the interpretation of these functions in the thermal reference states 
ujb = J dp{(3) ujp, there appears no problem with the energy density which can be 
determined by the local observable 6^^{x). Hence its expectation values are 



ujB{E^n = jdp{P)E^^{P). (A.6) 
It is less obvious, however, that the expectation values of S*^, 

u:b{S^)= fdp{(3)S^{(3), (A.7) 



can be interpreted as entropy currents of the non-equilibrium states ub- In or- 
der to justify this interpretation, we make the following physically meaningful 
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assumptions: the entropy density 5 is a concave function of the energy density 
E = Q — P, and the pressure P as well as E — P increase monotonically with 
the temperature T = The latter condition means that the negative free 

energy (the isothermal work which can be extracted from equilibrium states in 
their rest systems) increases with increasing temperature, yet less rapidly than the 
total energy. 

We will show that these assumptions entail a version of the second law based 
on ( I A .71 ) . To this end let us determine in any given Lorentz frame e' G K|. , e' ^ = 1 , 



the states lob G C with fixed energy density e'^e^ujBiE^'^) = E' which maximize 
the entropy functional ojb ^ ^'ij!^b{,S^). Considering first the KMS states and 
making the substitution (3 = (/3^)^/^e — > iE,e), we have 

e'/u = {{E + Pm {e'ef - P{E)) = E' (A.8) 

e'^u;p{S'^) = S{E){e'e). (A.9) 



Inserting (e'e) from relation ( |A.8| ) into ( [A. 91 ), we get 

e'^uj.iSn = S{E){{E' + P{E))/{E + P{E))Y^'. (A.IO) 

The logarithmic derivative of the right hand side of this equality with respect to 
E satisfies 

{dS/dE)/S +{dP/dE)/2{E' + P) - (1 + dP/dE)/2{E + P) > 

1/ST- {l + dP/dE)/2{E + P) = {l - dP/dE) /2{E + P) > 0, ^^''^^^ 



where we made use of (|A.3| ) and of the preceding assumptions, implying < 



dP/dE < 1. Thus the expression ( |A.10D attains its maximum for the maximal 
possible energy E = E' (corresponding to e = e') compatible with (|A.8|), i.e. 



e'^oop{Sn = S{E){e'e)<S{E'). (A.12) 
Turning to the general case, let ub = f dp{(3) ojp be any state with 

e'/,ojB{E^n = I da{E, e) {{E + P{E)) {e'ef - P{E)) = E' , (A.13) 

where a is the normalized measure obtained from p by the substitution j3 — > {E, e). 
In view of the preceding results and the concavity of S we obtain 

e'^LUBiS") = j da{E,e)SiE) (e'e) < jda{E,e)S{iE + P{E)) {e'ef -P{E)) 

< s[jda{E,e){{E + P{E)){e'ef-P{E)))=S{E'), (A.14) 



and equality holds in ([A.14|) if and only if the measure a is concentrated at the 



point E',e'. Thus ujpi, where (3' = (/5'^)^/^e' and P'"^ is fixed by E', is the unique 
state in C which maximizes the entropy functional ujb ^ ^'^i^b{S^) under the 

27 



given conditions. As this functional is also additive, the interpretation of ( |A.7| ) as 
entropy current is justified. 

Similarly, one can define the mean free energy density of the states in C, 

UB{Fn = I dpiP) = - [ dpi(3) P{f3') g^r (A.15) 



It thus coincides with the negative mean pressure in all Lorentz frames and has a 
sharp value for all states ujb with a definite temperature T, corresponding to mea- 
sures p which are concentrated on the manifold (3"^ = T~^. Hence the set of states 
minimizing the free energy functional ujb ^ ^'^^'u^b{F^^) for given temperature 
is degenerate. 

We conclude this appendix by illustrating the preceding notions for the case 
of a dilation invariant theory, where the stress energy tensor 9^"{x) is trace- 
less, (x) = 0, and transforms canonically under the dilations according to 
6s{9^'^{x)) = s^O^"{sx), s G M_|_. A special example is the model discussed in 
Sec. 5. It follows from the definition of KMS states and their uniqueness that 
uJi3 o Ss = cUs-i/j. Plugging these pieces of information into equation ( |A.1| ), we 
obtain 

E'^'iP) = C (Ap^p" - P^g^"'){p^)-\ (A.16) 

where C is a constant depending on the underlying theory. Thus the temperature 
dependence of the energy density in the rest system is given by E{P) = 3 C(/5^)~^, 
in accordance with Stefan-Boltzmann's law, and the pressure satisfies 

P(/5) = C(/52)-2. (A.17) 

So E = 3P, which also follows directly from the fact that the stress energy tensor 
is traceless. In particular, < dP/dE < 1, as was anticipated in the preceding 
analysis. Finally, the entropy current is given by 

S^'^p) =AC p^'ip^)-^, (A.18) 

so the pertinent thermal functions are completely fixed apart from the value of the 
constant C . 
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